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Abstract 

Formulas are derived for solutions of many-body wave scattering 
problems by small particles in the case of acoustically soft, hard, and 
impedance particles embedded in an inhomogeneous medium. The lim- 
iting case is considered, when the size a of small particles tends to zero 
while their number tends to infinity at a suitable rate. Equations for 
the limiting effective (self-consistent) field in the medium are derived. 
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1 Introduction 

There is a large literature on wave scattering by small bodies, starting from 
Rayleigh's work (1871). For the problem of wave scattering by one body 
an analytical solution was found only for the bodies of special shapes, for 
example, for balls and ellipsoids. If the scatterer is small then the scattered 
field can be calculated analytically for bodies of arbitrary shapes, see [7], 
where this theory is presented. 

The many-body wave scattering problem was discussed in the literature 
mostly numerically, if the number of scatterers is small, or under the as- 
sumption that the influence of the waves, scattered by other particles on 
a particular particle is negligible ( see [5], where one finds a large bibli- 
ography, 1386 entries). This corresponds to the case when the distance d 
between neighboring particles is much larger than the wavelength A, and the 
characteristic size a of a small body (particle) is much smaller than A. By 
k = 3j£ the wave number is denoted. 

The basic results of our paper consist of: 
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i) Derivation of formulas for the scattering amplitude for the wave scat- 
tering problem by one small (ka <C 1) body of an arbitrary shape under the 
Dirichlet, impedance, or Neumann boundary condition (acoustically soft, 
impedance, or hard particle), 

ii) Solution to many-body wave scattering problem by many such par- 
ticles under the assumptions a <C d and a « A, where d is the minimal 
distance between neighboring particles, 

hi) Derivation of the equations for the limiting effective (self-consistent) 
held in the medium when a — > and the number M = M(a) of the small 
particles tends to infinity at an appropriate rate, 

iv) Derivation of linear algebraic systems for solving many-body wave 
scattering problems; these system are not obtained by a discretization of 
boundary integral equations. 

Let us formulate the wave scattering problems we deal with. First, let 
us consider a one-body scattering problem. Let D\ be a bounded domain in 
]R 3 with a sufficiently smooth boundary S\. The scattering problem consists 
of finding the solution to the problem: 

(V 2 + k 2 )u = 0m D' 1 :=R 3 \D 1 , (1) 

Pu = 0onSi, (2) 

u = u + v, (3) 

where 

Uo = e ika - x , a£S 2 , (4) 

S 2 is the unit sphere in M. 3 , uq is the incident field, v is the scattered field 
satisfying the radiation condition 

v r — iKv = o - , r := \x — >■ oo, v r := — , (5) 
\r J ' or 

Tu is the boundary condition (be) of one of the following types 

Tu = Tiu = u (Dirichlet be), (6) 

Tu = T2U = un — Ciu, LmCi ^ 0, ( impedance be), (7) 
where £1 is a constant, N is the unit normal to Si, pointing out of Di, and 

Tu = T3U = un, ( Neumann be). (8) 



2 



It is well known (see, e.g., [6]) that problem (P)-© has a unique solution. 
We now assume that 

a := 0.5diamDi, ka < 1, (9) 

and look for the solution to problem (HJ-Q of the form 

r e ik\x-y\ 

u(x) = u (x)+ g(x,t)ai(t)dt, g(x,y) := — f , (10) 

J si 47r k - y\ 

where dt is the element of the surface area of S\. One can prove that the 
unique solution to the scattering problem ([l])-® with any of the boundary 
conditions ©-([8]) can be found in the form (jlQj) . and the function u\ in 
equation (|10p is uniquely defined from the boundary condition ([2]). The 
scattering amplitude A(/3, a) = A(f3, a, k) is defined by the formula 

e tkr , s /1\ „ x 



A(P, a, k) + o - , t — y oo, /3 := — . (11) 
r \r J r 

The equations for finding a% are: 

g(s,t)a 1 (t)dt = -u (s), (12) 

Si 

uoN-duo^ — — --Ci/ g{s,t)a 1 (t)dt = 0, (13) 

Ao\ — CTi , jN 

u 0N + 2 = 0, (14) 
respectively, for conditions ©-([H]). The operator A is defined as follows: 

f d 

Equations (|12|) -(|14|) are uniquely solvable, but there are no analytic formulas 
for their solutions for bodies of arbitrary shapes. However, if the body D\ 
is small, ka -C 1, one can rewrite fllOf) as 



u(x) = u (x)+g(x,0)Qi+ \g(x,t) - g(x,0)]ai(t)dt, (16) 

J Si 

where 

Q x := [ a x {t)dt, (17) 

iSi 
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and € D\ is the origin. 

If ka <C 1, then we prove that 



\g(x,0)Q 1 \ » 



fl'(M) - ff(x,0)]cri(t)di 
Si 



\x\ > a. (18) 



Therefore, the scattered field is determined outside D\ by a single number 
Q\. This number can be obtained analytically without solving equations 
(fT2"j) - (fT3j) . The case (fT4"j) requires a special approach by the reason discussed 
in detail later. 

Let us give the results for equations (fT2"j) and (Q2J) first. For equation 
(fl~2j) one has 

Qi= [ a l (t)dt = -Cu Q (0) [l + o(l)], a^O, (19) 

where C is the electric capacitance of a perfect conductor with the shape 
D\. For equation f)13|) one has 



Qi = -C|£i|uo(0)[l + o(l)], a^O, (20) 

where |Si| is the surface area of S\. The scattering amplitude for problem 
(HI)-© with T = Ti (acoustically soft particle) is 

A 1 03,q) = -^[1 + O (1)], (21) 

since 

«o(0) = e ika -% =0 = 1. 

Therefore, in this case the scattering is isotropic and of the order 0(a), 
because the capacitance C = 0(a). 

The scattering amplitude for problem ([T])-© with r = T2 (small impedance 
particles) is : 

A 2 (a,/3) = -^[l + (l)], (22) 

since uq(0) = 1. 

In this case the scattering is also isotropic, and of the order 0((\Si\). 

If £ = 0(1), then A 2 = 0(a 2 ), because \Si\ = 0(a 2 ). If Ci = (±), 
k G (0, 1), then A 2 = 0{a 2 - K ). The case k = 1 was considered in [9]. 

The scattering amplitude for problem ([I])-© with r = (acoustically 
hard particles) is 

4,(0, a) = -^^(l + /3 pq P p a q ), iiu Q = e ika ' x . (23) 
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Here and below summation is understood over the repeated indices, a q = 
e q denotes the dot product of two vectors in M 3 , p, q = 1, 2, 3, {e p } 
is an orthonormal Cartesian basis of M 3 , \D\\ is the volume of D\, f3 pq is the 
magnetic polarizability tensor defined as follows ([7], p. 62): 

Ppi := tttt / tp"iq(f)dt> ( 24 ) 



\Di\ JSl 

o\ q is the solution to the equation 

a lq {s) = A a lq - 2N q (s), (25) 

N q (s) = N(s) ■ e q , N = N(s) is the unit outer normal to S\ at the point s, 
i.e., the normal pointing out of D\, and is the operator A at k = 0. For 
small bodies \\A — Aq\\ = o(ka). 

If uq(x) is an arbitrary field satisfying equation ([T]), not necessarily the 
plane wave e lka ' x , then 

A 3 a) = H£i (ikfa^Pp + Am ) • (26) 

The above formulas are derived in Section 2. In Section 3 we develop a 
theory for many-body wave scattering problem and derive the equations for 
effective field in the medium, in which many small particles are embedded, 
as a — > 0. 

The results, presented in this paper, are based on the earlier works of 
the author ([9]- [23]). Our presentation and some of the results are novel. 
These results and methods of their derivation differ much from those in the 
homogenization theory ([2], [1]). The differences are: 

i) no periodic structure in the problems is assumed, 

ii) the operators in our problems are non-selfadjoint and have continuous 
spectrum, 

iii) the limiting medium is not homogeneous and its parameters are not 
periodic, 

iv) the technique for passing to the limit is different from one used in 
homogenization theory. 

2 Derivation of the formulas for one-body wave 
scattering problems 

Let us recall the known result (see e.g., 0) 

° f g(x,t)a 1 (t)dt= Aai ~ ai (27) 
Si 1 



ON. 
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concerning the limiting value of the normal derivative of single-layer po- 
tential from outside. Let x m G D m , t 6 S m , S m is the surface of D m , 
a = 0.5 diam£) m . 

In this Section m = 1, and x m = is the origin. 

We assume that ka <C 1, ad^ 1 <C 1, so |cc — x m | = d S> a. Then 

ik(x-X m )°-(t-X m ) (l + + ) ) , (28) 



Air\x — 1\ &7r\x — x m \ V d 

k\x — t\ = k\x — x m \ — k(x — x m )° ■ (t — x m ) + O I — — ) , (29) 



d 



where 

j _ I _ | ( _ \o x ~ x m 

u — \x Xm | , yx x m ) . — . , 

and 



J£z!l =i + o(S). (30) 

Let us derive estimate ([19]). Since |f| < a on Si, one has 

s(s,t) =g (s,t)(l + O(ka)), 

where go(s,t) = 47r |*_ t | ■ Since uq(s) is a smooth function, one has |uo(s) — 
uo(0)| = 0(a). Consequently, equation (fT2|) can be considered as an equation 
for electrostatic charge distribution <J\(t) on the surface S\ of a perfect 
conductor D\, charged to the constant potential — Uq(Q) (up to a small term 
of the order 0(ka)). It is known that the total charge Q± = f s <j\(t)dt of 
this conductor is equal to 

Qi = -Cuo(0)(l + O(fca)), (31) 

where C is the electric capacitance of the perfect conductor with the shape 
D x . 

Analytic formulas for electric capacitance C of a perfect conductor of 
an arbitrary shape, which allow to calculate C with a desired accuracy, are 
derived in [7j. For example, the zeroth approximation formula is 

cm = -r™ ™ = ~ 8 '' (32) 

J Si JSi r st 



and we assume in (|32|) that eo = 1, where eo is the dielectric constant of the 
homogeneous medium in which the perfect conductor is placed. Formula 
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AMD is formula JT^. If u (x) = e** 0- *, then u (0) = 1, and Q 1 = -C(l + 
0(ka)). In this case 

A 1 ( / 9,a) = ^ = —[1 + 0(^)1, 

which is formula (|21|) . 

Consider now wave scattering by an impedance particle. 

Let us derive formula ([20]) . Integrate equation (Q2J) over Si, use the 
divergence formula 

/ u 0N ds = / V 2 u dx = -k 2 / u dx = £; 2 |L>i|u (0)[l + o(l)], (33) 
where \D\ \ = 0(a 3 ), and the formula 

-Ci / u ds = -Ci|fi ! i|uo(0)[l + o(l)]. (34) 
Futhermore \J S g(s, t)ds\ = 0(a), so 

Ci / ds g(s,t)ax(t)dt = 0(aQ x ). (35) 

JS! JSl 

Therefore, the term (|35p is negligible compared with Q\ as a — > 0. Finally, 
if fca <C 1, then g(s, t) = go(s, t) (1 + ik\s — t\ + . . .) , and 

g(s,t) = —g (s,t)[l + O(ka)]. (36) 



Denote by Aq the operator 

^ = 2/ ^(t)& (37) 
It is known from the potential theory that 

/ A a 1 ds = - f ai(t)dt, 2 I &m }-^ t ' ds = -1, teSi. (38) 

>/ S~\_ J S\ J Si s 

Therefore, 

f ds Aai ~ ai =-Q 1 [l + 0(ka)]. (39) 
J Si 2 

Consequently, from fromulas (|33l) - (|39l) one gets formula (|22|) . 
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One can see that the wave scattering by an impedance particle is isotropic, 
and the scattered field is of the order 0(Ci\Si\). Since |£i| = 0(a 2 ), one 
would have 0(&|Si|) = 0{a 2 ' K ) if Ci = O (£), « G (0, 1). 

Consider now wave scattering by an acoustically hard small particle, i.e., 
the problem with the Neumann boundary condition. 

In this case we will prove that: 

i) The scattering is anisotropic, 

ii) It is defined not by a single number, as in the previous two cases, but 
by a tensor, 

and 

iii) The order of the scattered field is 0(a 3 ) as a — > 0, for a fixed k > 0, 
i.e., the scattered field is much smaller than in the previous two cases. 

When one integrates over Si equation (fl~3l) . one gets 

(40) 



Qi= V 2 u dx = V 2 u (0)|£>i|[l + o(l)], a-)-0. 

JDi 



Thus, Qi = 0(a 3 ). Therefore, the contribution of the term e lkx0 ' 1 in for- 
mula ([5SD with 

%m — will be also of the order O(a^) and should be taken 
into account, in contrast to the previous two cases. Namely, 

u(x) = u (x) + g(x,0) [ e'^aU^dt, (3:=^- = x°. (41) 
J Si M 

One has 

/ e- lk ^ t a 1 (t)dt = Q 1 -ik(3 p [ t p o- x {t)dt, (42) 
JSi JSi 

where the terms of higher order of smallness are neglected and summation 
over index p is understood. The function o\ solves equation fjl4[) : 



o"i = Ao\ + 2uon = Aa\ + 2ika q N q uo(s), s € Si (43) 
if uq(x) = e tka ' x . 

Comparing ()43j) with (|25|) , using (|24p , and taking into account that ka <C 
1, one gets 



-ik(3 p / t p ai{t)dt = -ik(3p\Di\p pq (-ika q )uo(0)[l + 0{ka)] 

JSi (44) 

= -fc 2 |Di|/3 pg/ p a g « o (O)[l + 0{ka)}. 
From ([iO]). (|32|) and (jUJ) one gets formula p5|l, because V 2 n = -k 2 u . 
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If Uq(x) is an arbitrary function, satisfying equation (pQ), then ika q in 
(|4"5|) is replaced by and —k 2 uo = Auo, which yields formula (|26l) . 

This completes the derivation of the formulas for the solution of scalar 
wave scattering problem by one small body on the boundary of which the 
Dirichlet, or the impedance, or the Neumann boundary condition is imposed. 



3 Many-body scattering problem 

In this Section we assume that there are M = M(a) small bodies (particles) 
D m , 1 < m < M, a = 0.5maxdiamZ) m , ka <C 1. The distance d = d(a) 
between neighboring bodies is much larger than a, d 3> a, but we do not 
assume that d ^> A, so there may be many small particles on the distances 
of the order of the wavelength A. This means that our medium with the 
embedded particles is not necessarily diluted. 

We assume that the small bodies are embedded in an arbitrary large 
but finite domain D, D C M 3 , so D m C D. Denote D' := M 3 \ D and 
f2 := U^l =1 D m , S m := dD m , dQ = U^ =1 S m . By N we denote a unit normal 
to dfl, pointing out of O, by \D m \ the volume of the body D m is denoted. 

The scattering problem consists of finding the solution to the following 
problem 

(V 2 + k 2 )u = in M 3 \ ft, (45) 

Tu = on dn, (46) 

u = u + v, (47) 

where uq is the incident field, satisfying equation (|45]) in R 3 , for example, 
uo = e tka ' x , a G S 2 , and v is the scattered field, satisfying the radiation 
condition ©. The boundary condition (f4l)j) can be of the types ©-([8]). 
In the case of impedance boundary condition ([7]) we assume that 

un = CmU on S m , 1 < m < M, (48) 

so the impedance may vary from one particle to another. We assume that 

Cm = %2, k £ (0, 1), (49) 
a K 

where x m G D m is a point in D m , and h(x), x G D, is a given function, 
which we can choose as we wish, subject to the condition Imh(x) < 0. For 
simplicity we assume that h{x) is a continuous function. 

Let us make the following assumption about the distribution of small 
particles: if A C D is an arbitrary open subset of D, then the number 
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A/"(A) of small particles in A, assuming the impedance boundary condition, 
is: 

A/- C (A) = ^L_ f N(x)dx[l + o(l)\, a^O, (50) 
a J A 

where N(x) > is a given function. If the Dirichlet boundary condition is 
assumed, then 

A/"d(A) = — f N(x)dx[l + o(l)], a^O. (51) 
a J a 

The case of the Neumann boundary condition will be considered later. 

We look for the solution to problem (|45p - ([47p with the Dirichlet bound- 
ary condition of the form 



M 

u 



uo+J2 g(x,t)a m (t)dt, (52) 

m=l J Sm 

where cr m (t) are some functions to be determined from the boundary condi- 
tion (|4"Uj) . It is proved in [U] that problem (|4l)j) - (|4"?]) has a unique solution of 
the form (|52p . For any a m (t) function (|52p solves equation (|45p and satisfies 
condition (|47p . The boundary condition (I46p determines o~ m uniquely. How- 
ever, if M > 1, then numerical solution of the system of integral equations 
for o~ m , 1 < m < M, which one gets from the boundary condition (I46h . is 
practically not feasible. 

To avoid this principal difficulty we prove that the solution to scattering 
problem (|45p - (|47|) is determined by M numbers 

Qm ■= / cr m (t)dt, (53) 

rather than M functions o~ m (t). 

This is possible to prove if the particles D m are small. We derive ana- 
lytical formulas for Q m as a — > 0. 

Let us define the effective (self-consistent) field u e (x) = Ue\x), acting 
on the j— th particle, by the formula 

u e (x) := u(x) — / g(x,t)aj(t)dt, \x — Xj\ ~ a. (54) 



Physically this field acts on the j— th particle and is a sum of the incident 
field and the fields acting from all other particles: 

u e (x) = u{ j) (x) := u (x) + ^2 g(x,t)a m (t)dt. (55) 



10 



Let us rewrite (l55j) as follows: 

M M 



M M .. 

u e (x) = u (x) + ^2 g(x,x m )Q m + ^2 / [g(x,t)-g(x,x m )]a m (t)dt. (56) 

We want to prove that the last sum is negligible compared with the first one 
as a — > 0. To prove this, let us give some estimates. One has \t — x m \ < a, 

d — | x Xffi | , 

\g(x,t) -g(x,x m )\ = max jo (^j ,0 (~^) \ > \d(x,x m )\ = 0(l/d). 

(57) 

Therefore, if \x — Xj\ = 0(a), then 

f Sm [g{x,t) - g(x,x m )]a m (t)dt ^ _ 



< 0(ad~ L + ka). (58) 



| g(x, x m )Q rn | 
One can also prove that 

J 1 /J 2 = 0(ka + ad- 1 ), (59) 

where J\ is the first sum in (j56[) and J2 is the second sum in (|56|) . Therefore, 
at any point 1 £ fi' = I 3 \ (1 one has 

M 

u e {x) = u Q (x) + ^2 g(x,x m )Q m , xett', (60) 

m=l 

where the terms of higher order of smallness are omitted. 
3.1 The case of acoustically soft particles 

If (|46p is the Dirichlet condition, then, as we have proved in Section 2 (see 
formula (j3~T1) ). one has 

Qm — Cra"e(j'm)' (^1) 

Thus, 

M 

u e (x) = u Q (x) - ^2 g(x,x m )C m u e (x m ), x £ O'. (62) 

m=l 

One has 

u(x) = u e (x) + o(l), a 0, (63) 
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so the full field and effective field are practically the same. 

Let us write a linear algebraic system (LAS) for finding unknown quan- 
tities u e (x m ): 

M 

u e (xj) = u (xj) - ^2 g{xj,x m )C m u e (x m ). (64) 

rrijtj 

If M is not very large, say M = O(10 3 ), then LAS (}64"|) can be solved 
numerically, and formula ([62]) can be used for calculation of u e (x). 

Consider the limiting case, when a — > 0. One can rewrite fl64f) as follows: 



Ue(£q) =U (£,q) -^2g(£q,£p)u e (Q ^ C m , (65) 

where {A p }^ =1 is a union of cubes which forms a covering of D, 

m&xdiamAp := b = b(a) a, 

p 

lim 6(a) = 0. (66) 

a— 5>0 

By |A p | we denote the volume (measure) of A p , and £ p is the center of A p , 
or a point x p in an arbitrary small body D p , located in A p . Let us assume 
that there exists the limit 

\im Exm l A f m =C(Q, ? P GA, (67) 

a^U \^p\ 

For example, one may have 

C m = c(£ p )a (68) 

for all m such that x m € A p , where c{x) is some function in D. If all D m 
are balls of radius a, then c(x) = 4tt. We have 

J2 C m = C p aM(A p ) = C p N(C p )\A p \[l + o(l)}, a -+ 0, (69) 



so limit (|67|) exists, and 

c(e P ) = c (e P )ivfe). (70) 
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From ([BS]), ([55 ]) -([7D ]) one gets 



u e {Q = u (C q ) -Y,9(C q ^ P )c(QN(Qu e (Q\A p \, 1< P <P. (71) 

Linear algebraic system (|7ip can be considered as the collocation method for 
solving integral equation 

u(x) = u (x)- g(x,y)c(y)N(y)u(y)dy. (72) 
Jd 

It is proved in [23] that system (|71[) is uniquely solvable for all sufficiently 
small b(a), and the function 

P 

u p( x ) : = ^Xpi^UeUp) (73) 

p=l 

converges in L°°(D) to the unique solution of equation (|72p . The function 
Xp(x) in (|73p is the characteristic function of the cube A p : it is equal to 1 in 
A p and vanishes outside A p . Thus, if a — > 0, the solution to the many-body 
wave scattering problem in the case of the Dirichlet boundary condition is 
well approximated by the unique solution of the integral equation (|72p . 

Applying the operator Lq := V 2 + k 2 to (|72j) . and using the formula 
Log(x,y) = —5(x — y), where 8(x) is the delta-function, one gets 

V 2 u + k 2 u - q(x)u = in E 3 , q(x) := c(x)N(x). (74) 

The physical conclusion is: 

If one embeds M(a) = 0(1 /a) small acoustically soft particles, which are 
distributed as in ()51[) . then one creats, as a — >• 0, a limiting medium, which 
is inhomogeneous, and has a refraction coefficient n 2 (x) = 1 — k~ 2 q(x). 

It is interesting from the physical point of view to note that the limit, 
as a — > 0, of the total volume of the embedded particles is zero. 

Indeed, the volume of one particle is 0(a 3 ), the total number M of the 
embedded particles is 0(a 3 M) = 0(a 2 ), and lim a _>.o G(a 2 ) = 0. 

The second observation is: if (|5ip holds, then on a unit length straight 
line there are 0(^73) particles, so the distance between neighboring particles 
is d = 0(a 1//3 ). If d = 0(a 7 ) with 7 > |, then the number of the embedded 
particles in a subdomain A p is O(jp) = 0(a~ 37 ). In this case, for 37 > 1, 
the limit in (|69j) is C(£ p ) = lim^o c p aO(a~^) = 00. Therefore, the product 
of this limit by u remains finite only if u = in D. Physically this means 
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that if the distances between neighboring perfectly soft particles are smaller 
than 0(a 1 ^ 3 ), namely, they are 0(a 7 ) with any 7 > |, then u = in D. 

On the other hand, if 7 < |, then the limit C(£ p ) = 0, and u = uq in D, 
so that the embedded particles do not change, in the limit a — > 0, properties 
of the medium. 

This concludes our discussion of the scattering problem for many acous- 
tically soft particles. 



3.2 Wave scattering by many impedance particles 



We assume now that (j49j) and (|50j) hold, use the exact boundary condition 
with r = T 2 , that is, 

UeN ~ CmUe + m<Tm — - Cm f ff( s > t)<J m {t)dt = 0, (75) 



and integrate (|75p over S m in order to derive an analytical asymptotic for- 
mula for Q m = J s a m (t)dt. 
We have 

/ u eN ds = i V 2 u e dx = 0(a 3 ), (76) 

JSm Jd 

■in 

Q m u e (s)ds = h(x m )a~ K \S m \u e (x m )[l + o(l)], a ->■ 0, (77) 

I' A m a m -a m ds = _ Qm[i + o{i)l a ^ Q) (78) 

J Sm. 



and 

Cm 



/ / g(s,t)a m {t)dt = h{x m )a 1 K Q m = o(Q m ), < k < 1. (79) 

Prom (|?5 > (f?9~ l) one finds 

Q m = -/i(x m )a 2 - K |5 m |a-V(^ m )[l + o(l)]. (80) 

This yields the formula for the approximate solution to the wave scattering 
problem for many impedance particles: 

M 

u(x) = u (x) - a 2 ~ K g(x, x m )b m h(x m )u e (x m )[l + o(l)], (81) 

m=l 

where 
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are some positive numbers which depend on the geometry of S m and are 
independent of a. For example, if all D m are balls of radius a, then b m = 4ir. 
A linear algebraic system for u e (x m ), analogous to (|64p . is 

M 

u e (xj) = u (xj) - a 2 ~ K ^2 g{xj 1 x m )b rn h{x m )u e (x m ). (82) 

The integral equation for the limiting effective field in the medium with 
embedded small particles, as a — > 0, is 

u(x) = u Q (x) -b g(x,y)N(y)h(y)u(y)dy, (83) 
Jd 

where 

u(x) = limit e (x), (84) 

a— »0 

and we have assumed in f)83[) for simplicity that b m = b for all m, that is, 
all small particles are of the same shape and size. 

Applying operator Lq = V 2 + k 2 to equation ([83]) . one finds the differ- 
ential equation for the limiting effective field u{x): 

(V 2 + k 2 - bN(x)h(x))u = in M 3 , (85) 

and u satisfies condition (j47j) . 

The conclusion is: the limiting medium is inhomogeneous, and its prop- 
erties are described by the function 

q(x) := bN(x)h(x). (86) 

Since the choice of the functions N(x) > and h(x), lmh(x) < 0, is at our 
disposal, we can create the medium with desired properties by embedding 
many small impedance particles, with suitable impedances, according to the 
distribution law (]50p with a suitable N(x). The function 

1 - k- 2 q(x) = n 2 (x) (87) 

is the refraction coefficient of the limiting medium. Given a desired refrac- 
tion coefficient n 2 (x), Imn 2 (x) > 0, one can find N(x) and h(x) so that 
(|87p holds, that is, one can create a material with a desired refraction coeffi- 
cient by embedding into a given material many small particles with suitable 
boundary impedances. 

This concludes our discussion of the wave scattering problem with many 
small impedance particles. 
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3.3 Wave scattering by many acoustically hard particles 

Consider now the case of acoustically hard particles, i.e., the case of Neu- 
mann boundary condition. The exact boundary integral equation for the 
function a m in this case is: 

UeN H ^ = ° - 

Arguing as in Section 2, see formulas ()40|) - ()44|) . one obtains 



M 

u e {x) = u (x)+'^2 9{x,x m ) 

m=l 



'"' \x-x m \ d(x) 



I -Dr. 

(89) 

Here we took into account that the unit vector ft in (I44p is now the vector 



l^_^ m | , and j3 p = \ x _ x "[ i where (x) p := x ■ e p is the p— th component of 

p}p=i- 



vector x in the Euclidean orthonormal basis {e p } 3 



There are three sets of unknowns in ([89]) : u e (x m ), 9u q^)^ ■> an d Au e (x m ), 
1 < m < M, 1 < q < 3. To obtain linear algebraic system for u e {x m ) and 
H%x) one sets x = * n ^9]) - takes the sum in ([89]) with m ^ j. This 
yields the first set of equations for finding these unknowns. Then one takes 
derivative of equation ([59"]) with respect to (x) g , sets x = Xj, and takes the 
sum in (|89p with j. This yields the second set of equations for finding 
these unknowns. Finally, one takes Laplacian of equation (|89l) . sets x = Xj, 
and takes the sum in (|89p with m ^ j. This yields the third set of linear 
algebraic equations for finding u e (x m ), ^q^™^ , and Au e (x m ). 

Passing to the limit a — > in equation (|89p. yields the equation for the 
limiting field 

u(x) = u (x)+ [ g(x,y) ( p{y)V 2 u{y) + ~ Vp . B pq (y)) dy, (90) 

Jd V oy q \x-y\ J 

where p(y) and B pq (y) are defined below, see formulas ([92]) and ([93]) . 

Let us derive equation (|90p . We start by transforming the sum in (|89p. 
Let {A;}^ be a covering of D by cubes A;, max/ diamA/ = 6 = b(a). We 
assume that 

6(a) > d » a, lim 6(a) = 0. 

a->0 

Thus, there are many small particles D m in A;. Let xi be a point in A/. 
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One has 

M 

^2 g{x,x m ) 



m=l 



\D r . 



X, ) 



1=1 



(91) 



Assume that the following limit exist: 

lim ?^ eA < 



l-Dr 



A; | 

(m). 



p(z/)> 



and 



J™ nri = Bpqiv), 

o->0, J/GA; I A/ 1 



<9ii e (y) <9u(y) 



(92) 
(93) 



lim u e (y) = lim "7^ = lim V 2 ?i e (y) = V 2 u(y). (94) 

a-»0 a-)-0 0{y)g Oy q a^O 

Then, the sum in (|91h converges to 



D 



du(y) x p - y p 
dy g \x - y\ jL '"' /1 



y) ( p(y)V 2 u(y) + ik ^l^ 1 ^ ^ BpJy) ) dy. (95) 



Consequently, ([89]) yields in the limit a — )■ equation ([90]) , Equation ([90]) 
cannot be reduced to a differential equation for u(x), because ([90]) is an 
integrodifferential equation whose integrand depends on x and y. 



4 Scattering by small particles embedded in an 
inhomogeneous medium 

Suppose that the operator V 2 + k 2 in ([1]) and in ([43]) is replaced by the 
operator Lq = V 2 + /c 2 ng(x), where ng(x) is a known function, 



Imng(x) > 0. 



(96) 



The function Uq(x) is the refraction coefficient of an inhomogeneous medium 
in which many small particles are embedded. The results, presented in Sec- 
tion 1-3 remain valid if one replaces function g(x, y) by the Green's function 
G(x,y), 

[V 2 + k 2 n 2 (x)]G(x, y) = -5(x - y), (97) 
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satisfying the radiation condition. We assume that 

njj(x) = 1 in D' := R 3 \ D. (98) 

The function G(x,y) is uniquely defined (see, e.g., [9]). The derivations of 
the results remain essentially the same because 

G(x,y) = g (x,y)[l + O(\x-y\)}, \x - y\ 0, (99) 

where (70(^1 y) = i 1T \ x-y \ • Estimates of G(:r, y) as |x — y| — >• and as |x — y| — )• 
00 are obtained in [S]. Smallness of particles in an inhomogeneous medium 
with refraction coefficient tlq(x) is described by the relation kn^a <C 1, where 
no := maxjigc |no(x)|, and a = maxi< m <MdiamZ) m . 
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